Charmonium suppression at RHIC: Signature of a strongly-interacting QGP, 

not a weakly interacting 
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Following a recent work on equation of state for strongly interacting quark-gluon plasma pj|, 
we revisited the equation of state by incorporating the nonperturbative effects in the deconfined 
plasma phase. Our results on thermodynamic observables viz. pressure, energy density, speed of 
sound etc. nicely fit with the lattice equation of state for gluon, massless and as well massive 
flavored plasma. Motivated by this agreement with lattice results, we have employed our equation 
of state to estimate the quarkonium suppression in an expanding, dissipative strongly interacting 
QGP produced in relativistic heavy-ion collisions and our prediction matches exactly with the recent 
PHENIX data on the centrality dependence of J/ip suppression in Au+Au collisions at BNL RHIC. 
We have also predicted for the T suppression in Pb+Pb collisions at LHC energy which could be 
tested in the ALICE experiments at CERN LHC. 
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I. INTRODUCTION 

Quantum chromodynamics (QCD) at high tem- 
perature is believed to be in quark gluon plasma 
(QGP) phase, whereby color charges are screened 
rather than confined. Asymptotic freedom of non- 
abelian gauge theories insures that for high enough 
temperature T 3> Aqcd, QGP is weakly coupled 
with dressed quarks and gluons behaving as quasi- 
particles near the ideal gas limit. 

In relativistic heavy ion collisions at RHIC, two 
novel phenomena - cone and ridge, notpresent in 
pp or d+Au collisions, were observed [J. Quark- 
gluon plasma which is produced in heavy ion colli- 
sions is not an ideal gas of quarks and gluons, it is 
rather a liquid having very low shear viscosity to en- 
tropy density (rj/s) ratio This strongly suggest 
that QGP may lie in the non-pcrturbative domain 
of QCD which is very hard to address both analyti- 
cally and computationally. Similar conclusion about 
QGP have been reached from recent lattice studies 
which predict that the equation of state (EoS) is 
interacting even at T ~ 4T C 0-Q3. Why QGP is 
strongly interacting in this temperature range and 
what is its meaning are not very well understood till 
the date. Since then, attempts have been made to 
understand strongly interacting nature of QGP and 
its small rj/s ratio, argued from large elliptic flow 
observed in RHIC data. Similar conclusions about 
the near perfect fluidity of QGP have been reached 
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from the AdS/CFT studies @, spectral functions 
and transport coefficients in lattice QCD [Tl| and 
studies based on classical strongly coupled plasmas 

There are several attempts to model the EoS of 
such strongly-interacting matter viz. bag model, 
confinement models, quasi-particle models, strongly 
interacting quark gluon plasmapl etc. In the 
bag model [l5f . QGP is treated as a big hadron with 
large number of partons interacting weakly but con- 
fined by the bag wall. Further inclusions of glue 
balls or hadrons improve the predictions. The con- 
finement models are the extension of bag model with 
smooth potential like Cornell potential [16| etc. with 
a better predictions. There is an interesting attempt 
[T7l [l8| to determine the EoS of interacting quarks 
and gluons upto Q(g 5 ) which has been further im- 
proved upon [191 . [20j by incorporating the contribu- 
tions from the nonperturbative scales viz. gT and 
g 2 T up to 0(.9 6 ln(l/ 5 ))[2H]. On the other hand, 
a semiclassical approach aimed to study the bulk 
properties of QGP automatically incorporate hard 
thermal loop (HTL) effects[22], [23| where the non- 
perturbative features manifested as effective mean 
color fields 0]. 

There are different versions of quasi-particle mod- 
els [25| where equation of state was derived with 
temperature dependent parton masses and bag 
constant [H, [27j . with effective degrees of free- 
dom [28|], etc. All of them claim to explain lat- 
tice results, either by adjusting free parameters in 
the model or by taking lattice data on one of the 
thermodynamic quantity as an input and predict- 
ing other quantities. However, physical picture of 
quasi-particle model and the origin of various tem- 
perature dependent quantities are not clear yet [29| . 
In strongly interacting QGP [30|-[32j], one considers 
all possible hadrons even at T > T c and try to ex- 
plain non-ideal behavior of QGP near T c . Recently, 
an equation of state for strongly-coupled plasma has 
been inferred by utilizing the understanding from 
strongly coupled QED plasma [33| which fits lattice 
data well. To be honest, deep and comprehensive 
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understanding is still missing and future investiga- 
tions (both theoretical and forthcoming experiments 
at LHC) may throw more light on this very exciting 
and complex issue. 

A suppressed yield of quarkonium in the dilcp- 
ton spectrum, measured in experiments [34l [35| was 
proposed as a signature of QGP formation. Under- 
standing the experimental data, however, turned out 
to be more complicated because the suppression pat- 
tern seen is not only due to the hot medium effects 
of screening [TTI . [36| , but more likely due to inter play 
of different interactions viz. cold nuclear matter [37| 
effect, gluonic dissociation [38| as well as those of 
recombination [39| of cc pairs etc. In order to disen- 
tangle these different effects we must know the prop- 
erties of quarkonium in medium and their dissocia- 
tion. There are two main lines of theoretical studies 
to determine quarkonium spectral functions at fi- 
nite temperature: potential models [4(|-|47j which 
have been widely used to study quarkonium, but 
their applicability at finite t emp erature is still un- 
der scrutiny and lattice QCD [UESl which provides 
the most straightforward way to determine spectral 
functions, but the results suffer from discretization 
effects and statistical errors, and thus are still incon- 
clusive. 

The short and intermediate distance properties 
of the heavy quark interaction are important for 
the understanding of in-medium modifications of the 
heavy quark bound states whereas the large distance 
behavior plays a crucial role in understanding the 
bulk properties of the QCD plasma phase [4f|. In the 
study of bulk thermodynamic quantities e.g. pres- 
sure, energy density etc., deviations from perturba- 
tive calculations and ideal gas behavior are found at 
temperatures much larger than the deconfinement 
temperature which is further supported by robust 
collective flows, strong jet and charm quenching, 
and charm flow, in RHIC data. This calls for a 
quantitative nonperturbativc calculations. Since the 
phase transition in full QCD appears as a crossover 
rather than a true phase transition [49], it is rea- 
sonable to assume that the string tension does not 
vanish abruptly above T c . So one should study its 
effects on the behavior of quarkonia in a hot QCD 
medium. Recently we [5(|l5l| had considered this 
potentially interesting issue by correcting the full 
Cornell potential with a dielectric function embody- 
ing the effects of the deconfined medium and not 
only its Coulomb part as usually done in the liter- 
ature. This led to a long-range Coulomb potential 
in addition to the usual Debye-screened form. With 
such an effective potential, we had investigated the 
effects of perturbative as well as nonperturbative ef- 
fects in QGP on the dissociation of different quarko- 
nium states H3, [5l| . 

Now let us consider a central collision in a nucleus- 
nucleus collision, which results in formation of QGP 
at initial time Tj. We assume the plasma to cool, 
according to Bjorken's boost invariant longitudinal 



expansion. There are many attempts [52M54| to in- 
corporate the plasma expansion dynamics to study 
charmonium suppression in relativistic nuclear col- 
lisions. However, some important points were not 
addressed in their works [52h54| to quantify the sup- 
pression in an expanding system: the viscous forces 
in the energy-momentum tensor and the proper EoS. 
The effects of dissipative forces were not included in 
hydrodynamic expansion which cause the plasma to 
expand slowly and results in more suppression. The 
equation of state employed in their works was ei- 
ther ideal or bag model EoS used to calculate two 
vital factors - screening energy density, e s of the 
system (corresponding to dissociation temperature) 
and the time, t s elapsed by the system to reach at 
e s , through expansion. Since the matter formed at 
RHIC is far from its ideal limit even at T > T c so the 
ideal or bag model equation of state is not reliable 
to study the suppression of charmonium yields in an 
expanding strongly interacting system. 

Recently we [55| had studied the survival of char- 
monium states in a dissipative strongly interacting 
QGP where the suppression of J/t/j's due to poten- 
tial screening in the deconfined medium at relativis- 
tic nuclear collisions is two step learning. First one is 
the understanding of dissociation in a static thermal 
medium for which knowledge of medium dependence 
of heavy quark potential is very much needed. Sec- 
ond one is the generalization of dissociation in an 
expanding medium where a equation of state plays 
the major role as an input. In [5a |. we used the 
equation of state for QGP, in analogy with strongly- 
coupled plasma where hadrons exist for T < T c and 
go to plasma of quarks and gluons (QGP) for T > T c 
and there is no hadrons or glue balls because confine- 
ment interactions due to QCD vacuum was assumed 
to be melted although Zahed and Shuryak [3l[ 
suggested that QGP at temperatures up to few T c 
supports weakly bound mesonic states. So the only 
interaction present in the deconfined plasma phase 
is Coulomb interaction and hence plasma parameter 
(r) becomes the ratio of average Coulomb poten- 
tial energy (=| ^-) to average kintic energy (~ T). 
Finally, an expression for equation of state was ob- 
tained as a function of plasma parameter, T 

Using the abovementioned equation of state, 
we [5|| got a better agreement with the PHENIX ex- 
perimental results [561 ] compared to earlier works |5~3j| 
but still it lacks complete agreement. The reason of 
disagreement may be two fold. This could be either 
due to the arbitrariness in the definition of dissoci- 
ation criteria or due to the equation of state which 
is not fully compatible with the nonperturbative na- 
ture of QGP. In the present work^we revisited the 
abovementioned equation of state [l| by incorporat- 
ing nonperturbative effects in the plasma parameter 
around which the equation of state is expanded. 

As discussed above, the existence of nonpertur- 
bative interactions even at T > T c indicates that 
the string tension may not vanish abruptly at T c , 
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so potential in the deconfincd phase could have a 
nonvanishing confining (string) term, in addition to 
the Coulomb term |50l | unlike Coulomb interaction 
alone in the aforesaid model [J. This is the central 
theme of our work where EoS has been obtained by 
retaining both terms in the potential and then cal- 
culate the thermodynamic variables viz pressure, en- 
ergy density, speed of sound etc. Our results match 
nicely with the lattice results of gluon pj, 2-flavor 
(massless) as well as 3- flavor (massless) QGP [57j ] . 
There is also an agreement with (2+1) (two mass- 
less and one is massive) and 4 flavoured lattice re- 
sults too. Motivated by the agreement with lattice 
results, we employ our equation of state to study 
the J/ip suppression in an expanding plasma in the 
presence of viscous forces with the universal ratio 
rj/s = 1/47T @. We have found complete agreement 
with the PHENIX results on the centrality depen- 
dence of J/tj) suppression at RHIC energy. Further 
we have predicted the T suppression which could 
be the potential candidate for the LHC experiment 
where J / tp suppression could be marred by enhance- 
ment due to recombination of abundant cc pairs. 

The paper is organized as follows. In Sec. II. A, 
we briefly discuss our recent work on medium modi- 
fied potential. Using this effective potential, we have 
then developed the equation of state for strongly in- 
teracting matter and have shown our results on pres- 
sure, energy density and speed of sound etc. along 
with the lattice data in Sec. II. B. In Sec. III. A, we 
have employed the aforesaid equation of state to 
study boost-invariant (1+1) dimensional longitudi- 
nal expansion in the presence of viscous forces and 
estimate the survival probability in a longitudinally 
expanding QGP in Sec. III. B. Results and discussion 
will be presented in Sec. IV and finally, we conclude 
in Sec.V. 



II. EQUATION OF STATE FOR 
STRONGLY-INTERACTING QGP 

In this section, first we briefly discuss the medium 
modified effective potential in Sec. II. A. In Sec. II. B, 
we revisit the EoS of strongly coupled QGP devel- 
oped by Bannur [T| and then obtain our EoS as 
a function of plasma parameter obtained from the 
medium-modified potential, discussed in Sec. II. A. 



A. Medium modified effective potential 

In thcrmodynamical studies of QCD plasma 
phase, deviations from perturbative calculations and 
ideal gas behavior at temperatures much larger than 
the deconfinement temperature calls for quantita- 
tive non-pcrturbative calculations [4||. In light of 
this finding, one can not simply ignore the effects 
of string tension between the quark- antiquark pairs 
beyond T c . Recently, this issue had successfully 



been addressed in the context of the dissociation of 
quarkonium in QGP [5£j [5l[ where we assumed the 
potential between a heavy quark-antiquark at T = 
as the Cornell potential : V(r) = ~a/r + ar and 
then corrected its Fourier transform (FT) V(k), to 
incorporate the medium modifications, as 



V(k) 



V(k) 
e(k) 



(1) 



where e(k) is the dielectric permittivity [58| . given 
by 
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and V(k) is the Fourier transform (FT) of the Cor- 
nell potential: 



V{k) = -^(2M-- 



4a 



(3) 



Substituting Eqs.([2]) and ([3]) into (JlJ and then eval- 
uating its inverse FT one obtains the independence 
of the medium modified potential 



V{r) 



2<T 

2a 



2a 
m D 
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(4) 



where constant terms arise from the basic computa- 
tions of real-time static potential [59| or from real- 
and imaginary-time correlators in a thermal QCD 
medium[6(| and are introduced to yield the cor- 
rect limit of V(r,T) as T -> 0. The potential gj) 
thus obtained has an additional long range Coulomb 
term with an (reduced) effective charge in addition 
to the conventional Yukawa term. It is worth to 
mention [6l"l ] that one-dimensional Fourier transform 
of the Cornell potential in the medium yields the 
screened form as used in the lattice QCD to study 
the quarkonium properties which assumes the one- 
dimensional color flux tube structure. However, at 
finite temperature that may not be the case since the 
flux tube structure may expand in more dimensions 
[Tl| . Therefore, it is better to consider the three- 
dimensional form of the medium modified Cornell 
potential. 

Recently we had employed the above medium- 
modified potential to estimate the dissociation pat- 
tern of the charmonium and bottomonium states 
and also explore how the pattern changes as w e g o 
from perturbative to nonperturbative regime [50| |. 
The results obtained in perturbative to nonperturba- 
tive domain were closer to the results obtained from 
the study of spectral function constructed in lattice 
QCD and potential model studies, respectively. This 
medium modified potential will be employed to de- 
velop the equation of state in next section. 
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B. Equation of state 

The equation of state for the quark matter pro- 
duced in relativistic nucleus-nucleus collisions is an 
important observable and the properties of the mat- 
ter are sensitive to it. The expansion of QGP is quite 
sensitive to EoS through the speed of sound which, 
in turn, explores the sensitivity of the quarkonium 
suppression to the equation of state [53|, [54| ■ 

Lattice results Q show that bulk thermodynamic 
observables such as pressure and energy density de- 
viate from ideal gas behavior even at ~ 4 T c and 
approach the ideal limit very slowly. This sug- 
gest that there are still (nonperturbative) interac- 
tions present in the deconfmed medium. There have 
been many attempts to explain this strongly inter- 
acting matter formed during ultrarelativistic colli- 
sions. Many models^such as bag model[l5|, quasi- 
particle model [62|, [63( , strongly interacting quark- 
gluon plasma (QGP) model |14| were used to ex- 
plain the complicated matter formed. None of the 
models were fully accepted by the physics commu- 
nity even though some models were successful to a 
great extent. An equation of state for such strongly 
interacting matter was developed by treating QGP 
near and above T c in Cornell potential as Coulomb 
plus linear confinement term and compared with the 
lattice results for pure gauge, two-flavor, and three- 
flavor QGP [53 . Arnold and Zhai [l7| and Zhai and 
Kastening [18f have developed a perturbative EoS 
of interacting quarks and gluons upto 0(g 5 ) which 
has further been extended by Kajantie et al. [i~9H2C| 
by including the nonperturbative contributions viz. 
O(gT) and 0{g 2 T) up to 0{g 6 ln(l/p)) (2l|. On the 
other hand, a semiclassical approach is devised to 
study the bulk properties of QGP [H H [MS] 
where nonperturbative effects manifested as effec- 
tive mean color fields. These color fields have the 
dual role of producing the soft and semisoft partons, 
apart from modulating their interactions. The emer- 
gence of such effective field degrees of freedom, to- 
gether with a classical transport has been indicated 
earlier [6~H |65|. 

Recently Bannur [l| developed an equation of 
state for a strongly-coupled QGP by appropriate 
modifications of strongly-coupled plasma in QED to 
take account color and flavor degrees of freedom with 
the running coupling constant and a reasonably good 
fit to the lattice results was obtained. Let us briefly 
discuss the strongly coupled plasma in QED where 
the equation of state is expressed as a function of 
plasma parameter T [67| : 



"QED 



+ u ex (T) nT 



where the first term represents the ideal contribution 
and the deviations from ideal EoS is given by, 
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where uf^ e is given by 
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(7) 



which was derived by Abe [68| in the formalism of 
giant cluster expansion with the Euler constant 7 
and is valid for T < .1. The term u® x CP is given by 

!i° CP = -0.898004r + 0.96786r 1/4 

+ 0.220703r- 1/4 - 0.86097 , (8) 

which was numerically obtained for one component 
plasma and is valid all V < 180 [69fl . In the above 
equations, the plasma parameter T is the ratio of av- 
erage potential energy to the average kinetic energy 
Let us now consider strongly-coupled plasma in 
QCD where it was assumed that hadron exists for 
T < T c and goes to QGP for T > T c . That is, 
for T > T c , it is the strongly interacting plasma 
of quarks and gluons and no hadrons or glue balls 
because it was assumed that confinement interac- 
tions due to QCD vacuum has been melted p| at 
T = T c . Hence the only interaction present in the 
deconfmed plasma phase is the Coulomb interaction 
and so the plasma parameter (r) was evaluated as 
the ratio of average Coulomb potential energy to av- 
erage kinetic energy. As discussed earlier in Sec. II. A, 
we will retain confinement interactions, in addition 
to the Coulomb interaction, through the linear term 
in the potential which manifests in the nonzero val- 
ues of string tension even at T > T c in the potential 
Finally, the equation of state has been obtained 
by using the potential (4) in the plasma parameter 
after inclusion of relativistic and quantum effects as: 



3 + u ex (T)j nT 



(9) 



where the form of u ex (T) remains the same as in 
(6). In terms of ideal contribution, the scaled-energy 
density is written as 



e(r) 



e 
£sb 



1 



1 



where e SB is given by, 



3a f T 4 



(10) 



(11) 
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where a/ = (16 + 21 rif /2)tt 2 /90 is a constant which 
depends on degrees of freedom of quarks and gluons. 
Here we have employed the QCD running coupling 
in MS scheme [70j |. in compatible with lattice simu- 
lation, up to two-loop level: 
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where b = (33 - 2n / )/(487r 2 ) and bi = (153 - 
19n/)/(3847r 4 ), p and A^jg are the scale parameter 

and the renormalization scale in MS scheme, respec- 
tively. For, the EoS to depend on the renormaliza- 
tion scale, the physical observables should be scale 
independent. We circumvent the problem by trad- 
ing off the dependence on A^jg to a dependence on 
the critical temperature T c . 

p exp(7 B + c) = A M g(T) 
A M s(T)exp( 7£ + c) = 4ttA t , (13) 

where c is a constant depending on colors and 
flavors: c = (n c — 4nj In 4) / (22n c — n/) and 
7^=0.5772156. There are several ambiguities, asso- 
ciated with the renormalization scale A^jg, the scale 
parameter p which occurs in the expression for the 
running coupling constant a s . This issue has been 
discussed well in literature and a popular way out 
is the BLM criterion due to Brodsky, Lepage and 
Mackenzie [7l| . In this criterion, A^jg is allowed to 
vary between ttT and AnT [72|. For our purposes, 
we choose the renormalization scale A^jg close to the 
central value 2ttT c [73| for ?z /=0 and nT c for both 
?i/ ~2 and nf=3 flavors. 

It is worth to mention here that if the factor 

In (2 In ) 

\— \ — » is much smaller than 1 (< 1) then 

In — 



■21, 
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the above expression reduces to the expression used 
in [l|, Eq.(10)], after neglecting the higher order 
terms of the above factor. However, this possibil- 
ity does not hold good for the temperature ranges 
used in the calculation and cause an error in coupling 
which finally makes the difference in the results be- 
tween our model and Bannur model [l|. Finally, we 
get the energy density e(T) from Eq.(10)and using 
the thermodynamic relation, 



e = T^-P 
dT 



we get the pressure as 



dTT 2 e{T{r)) /T 3 



(14) 



(15) 



T 



where Pq is the pressure at some reference temper- 
ature To and has been fixed with the values of pres- 
sure at critical temperature for a particular system 
- gluon plasma, 2-flavor plasma etc. Once we know 
the pressure P and energy density e, the speed of 
sound c 2 (= ^) can be evaluated. 

In Fig. 1, we have plotted the variation of pressure 
(P/T A ) with temperature (T/T c ) for pure gauge, 2- 
flavor and 3-flavor QGP along with lattice results. 
In Bannur model [l| , for each system, g c and At are 
adjusted to get a good fit to lattice results. However, 
in our calculation, there is no quantity to be fitted 



for predicting lattice results. We have fixed Pq from 
the lattice data at the critical temperature T c for 
each system, separately. 

Once pressure, P{T) is obtained, then other 
macroscopic quantities such as energy density e, 
speed of sound c 2 etc. can be derived from P{T) 
and no other parameters are needed. In Fig. 2, we 
plotted the energy density (s/T 4 ) with temperature 
(T/T c ) for all three systems along with lattice re- 
sults and a reasonably good fit is obtained without 
any extra parameters. All three curves looks similar, 
but shifts to left as flavor content increases. For the 
sake of comparison with the results of Bannur EoS, 
we have taken the critical temperatures T c equal to 
275, 175 and 155 MeV for gluon plasma, 2-flavor and 
3-flavor QGP, respectively. 

In Fig. 3, the speed of sound, c 2 is plotted for all 
three systems, but matching have been checked with 
the only available lattice results for gluon plasma. 
There is an excellent agreement with the lattice re- 
sults and we have predicted the results for the fla- 
vored QGP. All three curves have similar behaviour, 
i.e, sharp rise near T c and then flatten to the ideal 
value (1/3). However, c 2 is larger for larger flavor 
content of plasma. In the vicinity of critical tem- 
perature, fits or predictions may not be good, es- 
pecially for energy density e and c 2 which strongly 
depends on variations of pressure P with respect to 
temperature T. Lattice data also has large error 
bars very close to T c . However, except for small re- 
gion at T = T c , our results are very good for all 
regions of T > T c . It is interesting to note that 
recently Peshier and Cassing [74| also obtained sim- 
ilar results on the dependence of plasma parameter 
r in quasi-particle model and concluded that QGP 
behaves like a liquid, not weakly-interacting gas. 

Now we are interested to study for the realistic 
case where u and d quarks have very small masses (5- 
10 MeV) and strange quarks are having masses 150- 
200 MeV and charm quark with mass 1.5 GeV. Let 
gj counts the effective number of degrees of freedom 
of a massive Fermi gas. For a massless gas we have, 
of course, <// = n/. In general we define 



where, 



f=u,d,--- 



(16) 
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dx 
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In Fig. 4, we have shown our results on (2+1)- 
flavors and 4-flavors QGP along with lattice 
data [75], EH and replotted the variation of P(T)/T 4 
with temperature T/T c for all systems. Similar plots 
for energy density e(T)/T 4 with temperature T/T c 
for all systems is replotted in Fig. 5. We have also 
compared with the resuts from Bannur model. 



n f =0 n r = 2 n=3 




FIG. 1: Plots of P/T 4 as a function of T/T c from Bannur EoS, our EoS and lattice results for (a) pure gauge 
(extreme left figure), 2-flavor QGP (middle figure) and 3-flavor QGP (extreme right figure). In each figure, solid 
line represents the results obtained from Bannur EoS, dashed line represents the results from our EoS and diamond 
symbols reDresent lattice results. 




FIG. 2: Plots of e/T 4 as a function of T/T c . The notations are the same as in Fig.l. 





0.4 
0.35 

0.3 
0.25 
V 0.2 
0.15 

0.1 
0.05 



T~H~ 



I . I 



I 1 I 1 I 1 I 1 I 1 I 




Bannur EoS 
our EoS 



I i I ' I ' I i I i I 





0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 
T/T 



FIG. 3: Plots of c 2 a as a function of T/T c from our EoS and Bannur EoS for pure gauge, 2-flavor QGP and 3-flavor 
QGP. 




FIG. 4: variation of P/T 4 as a function of T/T c for a) two massless and one massive (2+1), b) and c) for 4-fiavour 
QGP for two different masses, m/T=0A and 0.2, respectively. The notations are the same as in Fig.l. 




This indicates that in the presence of a heav- 
ier quark the deviations of the pressure from the 
ideal gas value is larger than in the massless limit. 
This is in qualitative agreement with the observa- 
tions. In view of the agreement with the results of 
lattice equation of state, our EoS is a right choice 
for the strongly-interacting matter possibly formed 
at RHIC to calculate the thermodynamical quanti- 
ties viz. screening energy density (e s ), the speed of 
sound etc. to study the hydrodynamical expansion 
of plasma and finally, to estimate the suppression of 
J /tp in nuclear collisions. 



III. SUPPRESSION OF J /if) IN A 
LONGITUDINALLY EXPANDING PLASMA 



In Sec. III. A, we study hydrodynamic boost- 
invariant Bjorkcn expansion in (1 + 1) dimension 
with the EoS discussed in Sec.IIB as an input. In 
addition, we explore the effects of dissipative terms 
up to first-order in the stress-tensor. Then we turn 
our attention to derive the J/?p survival probability 
for an expanding QGP, in Sec. III. B. 



A. Longitudinal expansion in the presence of 
dissipative forces 

In the presence of viscous forces, the energy- 
momentum tensor is written as, 



where 77 is the co-efficient of the shear viscosity and 
(V M M y ) is the symmetrized velocity gradient. 

In (1+1) dimensional Bjorkcn expansion in the 
first-order dissipative hydrodynamics, only one com- 
ponent 7r T ' T ' of the viscous stress tensor is non-zero, 
hence the equation of motion reads, 



Or 



P 



ir] 
3^2 



(20) 



The first term in the RHS is the same as in the 
case of zeroth-order (non-viscous) hydrodynamics 
and the second term is the correction arising from 
constant 77/s which causes the system to expand 
slowly compared to the perfect fluid, rj = 0. 

The solution of equation of motion (19) is ob- 
tained as, 



e(r)r 



+ 



4q 
3f2' 



-(!+<£) 



<n)T- +—2 
on 



const , 



(21) 



where the constant a is ( 2 ) T?Tj and the symbols, 
f 2 and f 2 are given by (1 — c 2 )r 2 and (1 — c 2 )"rf, re- 
spectively. The first term accounts for the contribu- 
tions coming from the zeroth-order expansion (ideal 
fluid) and the second term is the first-order viscous 
corrections. The above equation (21) will finally be 
coupled with the dissociation of a quarkonium in a 
static thermal medium to estimate the quarkonium 
suppression in relativistic nucleus-nucleus collisions. 



T"" = (e+ P )u^u v + g^p + tt^, (18) 

where the stress-energy tensor, it 11 " up to first-order 
is given by 



(19) 



B. Survival probability 

We now have all the ingredients to write down the 
survival probability. Chu and Matsui [52| studied 
the transverse momentum dependence (jpr) of the 
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survival probability by choosing the speed of sound 
c 2 = 1/3 (ideal EoS) and the extreme value c 2 = 0. 
This work was further generalized by invoking the 
various parameters for Au-Au collisions at RHIC 
in [53| to include the effects of realistic EoS in an 
adhoc manner by simply choosing a lower value of 
c 2 = 1/5. Instead of taking arbitrary values of c 2 
we tabulated the values of c 2 in Tables II and IV- 
VI corresponding to the dissociation temperatures 
calculated from our EoS [55| for charmonium and 
bottomonium states, respectively. Moreover, in the 
light of recent experimental finding from RHIC, one 
cannot ignore the viscous effects while studying char- 
monium suppression. Here, we address these issues. 

Let us take an initial energy-density profile on a 
transverse plane : 



0(R T - r) (22) 



where r is the transverse co-ordinate and Rt is the 
transverse radius of the nucleus. One can define an 
average energy density (e,) as 



J 2ir r dre(r;Ti) 



so that 



(l + j3)(e i ) ;/3 = l. 



The average initial energy density (e;) 
given by the modified Bjorkcn formula: 



dEn 



A T n \ dy h 



(23) 
(24) 

(Zi will be 

(25) 



y h =o 



where At is the transverse overlap area of the 
colliding nuclei and (dE? / dyh) yh =o is the trans- 
verse energy deposited per unit rapidity of outgo- 
ing hadrons. Both depend on the number of par- 
ticipants N part (H3| and thus provide centrality de- 
pendent initial average energy density (e;) in the 
transverse plane (Table I). For this purpose, we 
have extracted the transverse overlap area At and 
the pseudo-rapidity distribution dEx/drfh \ m _ [HO] 
at various values of number of participants N part . 
These dEr/dr/h \ Vl =0 numbers are then multiplied 
by a Jacobian 1.25 to yield the rapidity distribution 
dET i 'dyh \ Vh= Q which will be further used to calcu- 
late the average initial energy density from Bjorken 
formula (25). The scaling factor £ = 5 has been 
introduced in order to obtain the desired values of 
initial energy densities [8l| for most central collision 
which are consistent with the predictions of the self- 
screened parton cascade model [82| and also with 
the requirements of hydrodynamic simulation [8ll | 
to fit the pseudo-rapidity distribution of charged 
particle multiplicity dN c h/drj for various centralities 



TABLE I: Kinematic characterization of Au+Au colli- 
sions at RHIC [51 



Nuclei a/sjvat V part 


<ei> 


Rt 


(GeV) 


(GeV/fm 3 ) 


(fin) 


22.0 


5.86 


3.45 


30.2 


7.92 


3.61 


40.2 


10.14 


3.79 


52.5 


12.76 


3.96 


66.7 


15.69 


4.16 


83.3 


18.58 


4.37 


Au+Au 200 103.0 


21.36 


4.61 


125.0 


24.38 


4.85 


151.0 


27.37 


5.12 


181.0 


30.52 


5.38 


215.0 


34.17 


5.64 


254.0 


37.39 


5.97 


300.0 


41.08 


6.31 


353.0 


45.09 


6.68 



observed in PHENIX experiments at RHIC energy. 
Later we will discuss the centrality dependence of 
initial energy densities at LHC energy in Sec. IV. 

The (screening) time, r s when the energy density 
of the system drops to the screening energy density 
e s is estimated from Eq.(f2"Tj) as 



i(r) 



4a 
3?F 



4a 
3f? 



(26) 



where e^(r) = e(r,;r) and f 2 is (1 — c 2 )t 2 . The 
critical radius r 8 , is seen to mark the boundary of 
the region where the quarkonium formation is sup- 
pressed, can be obtained by equating the duration 
of screening T s (r) to the formation time tp = 7T"f 
for the quarkonium in the plasma frame and is given 
by: 



r s = R T (1 - A) 1/2 ® (1 - A) 
where A is given by 



(27) 



A 



1 4a 
'T t 3r7 



t_F 



(1+c 



' 4o 
3% 



(28) 



with l 2 F = (1 — c 2 )t 2 F . The quark-pair will escape the 
screening region and form quarkonium if its position 
vector r and transverse momentum pt are such that 



|r + r F p T /A</| > r s . 



(29) 
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Thus, if <f> is the angle between the vectors r and pt, 
then the above condition reduces to 

cos > [(r 2 s - r 2 ) M - r| p 2 T /M] /[2rr F p T ] , 

(30) 

which leads to a range of values of (j) when the 
quarkonium would escape. Now we can write for 
the survival probability of the quarkonium: 



r dr 



-<£max 



-<^max 



>P(r,p T ) 



2 71 



rdr P(r 7 px) 



I 

(31) 



where </>max is the maximum positive angle (0 < 
4> < tt) allowed by Ea. ((5D|> : 



'max 



7r if y < — 1 

cos -1 if — 1 < j/ < 1 , 
if y > 1 



(32) 



where 



y = [(r 2 s - r 2 ) M - r 2 p 2 T /M] /[2rr F p T ] , (33) 

and P is the probability for the quark-pair produc- 
tion at (r, Pt), in a hard collision which may be 
factored out as 

P(r,p T ) = f(r)g(p T ), (34) 
where we take the profile function f(r) as 



f(r) oc 



1 



R 2 T 



0(R 7 



(35) 



with a = 1/2. 

Often experimental measurement of survival prob- 
ability at a given number of participants (N n ) or 
rapidity (y) is reported in terms of the p T -integrated 
yield ratio (nuclear modification factor) over the 

range p^ m < p T < p™ ax whose theoretical expres- 
sion would be 



max 



(S(P T )) 



T 

pinax 



(36) 



In nucleus- nucleus collisions, it is known that only 
about 60% of the observed J/ip originate directly in 
hard collisions while 30% of them come from the de- 
cay of Xc and 10% from the decay of tp'. Hence, the 
p T -integrated inclusive survival probability of J/ip 
in the QGP becomes [H[5|. 



, inci, , dir. , dir, , dir. 

(5 )= 0.6(5 ) +0.3(5 > +0.1(5 ) 



(37) 



IV. RESULTS AND DISCUSSIONS 

Before displaying the results, let us discuss the 
physical understanding of quarkonium suppression 
due to screening in the deconfined medium produced 
in relativistic nucleus-nucleus collisions. This in- 
volves a competition of various time-scales involved 
in an expanding plasma. First one is the screen- 
ing time, t s as the time available for the hot and 
dense expanding system during which J/ip's are sup- 
pressed. Second one is the formation time of J/ip 
in the plasma frame (tF = 7Tf) which depends on 
the transverse momentum by which the cc pairs was 
originally produced. Third one is the cooling rate 
which depends on the speed of sound, c 2 through 
the equation of state. The screening time not only 
depends upon the screening energy density, e s but 
also depends on the speed of sound through equa- 
tion of state. The value of g s is different for different 
charmonium states and is calculated from the equa- 
tion of state and hence varies from one EoS to other. 
If g s > ej, initial energy density, then there will be 
no suppression at all i.e., survival probability, S(p T ) 
is equal to 1. 

More precisely, the screening time depends upon 
i) the screening energy density and the difference be- 
tween a given initial energy density and screening 
energy density e s : the more will be the difference 
the more will be the suppression, ii) the speed of 
sound: for c 2 less than the ideal limit (1/3), the rate 
of cooling will be slower which, in turn, makes the 
screening time larger for a fixed difference in (gj- g s ) 
and leads to more suppression, and iii) the rj/s ra- 
tio: an additional handle to explore the equation of 
state by controlling the expansion of the plasma. If 
the ratio is non-zero then the cooling will be slower 
compared to rj/s=0, so the system will take longer 
time to reach g s resulting the higher value of screen- 
ing time and hence more suppression compared to 
rj/s = 0. With this physical understanding we ana- 
lyze our results, (S(p T )) as a function of the number 
of participants Np art in an expanding QGP. 

In our analysis, we have employed the dissocia- 
tion temperatures for the charmonium states (J/i/j, 
Xc etc.) computed from the lattice QCD correla- 
tor studies [85] in Table II. The corresponding val- 
ues of screening energy densities, e s and the speed 
of sound c 2 calculated in our EoS are also listed 
which will be used as inputs along with the kine- 
matic data in Table I, to calculate (S(p T )}. We have 
shown the variation of p T -integrated survival prob- 
ability (in the range allowed by invariant p T spec- 
trum of J j tjj in Phenix experiment [561 ] ) with Np art 
at mid-rapidity in Fig. 6. The experimental data (the 
nuclear- modification factor Raa) are shown by the 
squares with error bars whereas circles represent se- 
quential suppression. It is found that our calcula- 
tion matches completely with the experimental re- 
sults. To see the importance of confinement inter- 
actions in the deconfined phase, we have also calcu- 
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TABLE II: Formation time (fm), dissociation tempera- 
ture To O, the speed of sound <? a and the screening 
energy density e s (GeV/ fm 3 ) for charmonium states, 
calculated both in our and Bannur EoS, respectively. 



State 


TF 


T D 


c 2 s (our) 


c^(BAN) 


e s (our) 


e s (BAN) 


J/V 


0.89 


2.1 


0.308 


0.275 


29.33 


32.05 




1.50 


1.12 


0.255 


0.214 


01.94 


02.36 


Xc 


2.00 


1.16 


0.261 


0.220 


02.28 


0.220 



lated the p T -integrated survival probability (denoted 
by diamonds) using the equation of state in Bannur 
model [J where there are no confinement interac- 
tions in the deconfined phase and have not found 
the agreement, as seen in our case. This difference 
is due to the fact that the screening energy density 
calculated in our EoS is smaller than the value ob- 
tained from Bannur model (as seen in Table II). The 
smaller value of screening energy density e s causes 
an increase in the screening time and results in more 
suppression to match with the PHENIX results at 
RHIC. 

At RHIC energy, J/ip yields have been resulted 
from a balance between annihilation of J ftp's due 
to hard, thermal gluons [86|, [87} along with colour 
screening [52, Ha| and enhancement due to coales- 
cence of uncorrelated cc pairs [89l - [9l| which are pro- 
duced thermally at deconfined medium, at the phase 
boundary during (statistical) hadronization [92l . [93| . 
However, recent PHENIX data do not show a fully 
confirmed indication of J / ip enhancement except for 
the fact that (p\) of the data and shape of rapidity- 
dependent nuclear modification factor Raa{u) [5a | 
show some characteristics of coalescence production. 

On the other hand, at LHC energy, there will 
be an abundance of thermally produced cc pairs 
so that one cannot make any concrete conclusion 
about the possible formation of QGP whereas the 
number of bb pairs produced in hot deconfined 
medium will be meagre (because of large bottom 
quark-antiquark masses) so that the competition 
between the suppression due to screening and 
dissociation both and the enhancement due to 
recombination may be unlikely. Therefore, it will be 
interesting to estimate the T suppression at LHC 
energy because the initial energy density at ALICE 
experiments will be high enough to suppress the T's 
sequentially making the bottomonium suppression 
an unambiguous signature. 

Before finding the centrality (or impact parame- 
ter) dependence of T suppression at LHC energy, it 
is necessary to know the centrality dependence of 
average initial energy density (e;) in terms of the 
number of participants N part at LHC energy. Re- 
cently Eskola et al [12] computed the dependence 
of initial energy density, gluon, quark and antiquark 



TABLE III: Kinematic characterization of Pb+Pb colli- 
sions at LHC H2 



Nuclei ^s N n N palt 
(GeV) 


<£i> 

(GeV/fm 3 ) 


Rt 
(fm) 


22.0 


217.97 


3.45 


30.2 


236.12 


3.61 


40.2 


253.15 


3.79 


52.5 


269.11 


3.96 


66.7 


291.45 


4.16 


83.3 


315.55 


4.37 


Pb+Pb 5500 103.0 


341.05 


4.61 


125.0 


370.49 


4.85 


151.0 


400.42 


5.12 


181.0 


433.18 


5.38 


215.0 


463.66 


5.64 


254.0 


507.04 


5.97 


300.0 


550.88 


6.31 


353.0 


600.39 


6.68 



numbers produced in ultrarclativistic heavy ion col- 
lisions with beam energy in parton saturation model. 
This model is based on the argument that the effects 
of all momentum scales can be estimated by per- 
forming the computation at the saturation momen- 
tum scales. The main emphasis of the study was at 
LHC and RHIC energies, although it gives reason- 
able good results at SPS too. The dependence of 
atomic number and beam energies of initial energy 
density, number density, temperature etc. are given 
by m 

d = 0.103GcVfm- 3 A a504 (Vi) - 786 , (38) 
m = 0.370 fm- 3 A°- 383 (^) - 574 , (39) 
Ti = 0.111 GeV ^4 126 (v / s) ' 197 (40) 

Using this model, we have first calculated the cen- 
trality dependence of the initial energy densities at 
RHlC energy and found a good agreement with the 
known values (Table I) obtained from the centr ality 
dependence of observed particle multiplicities [80j j. 
Inspired by the success of the parton saturation 
model at RHIC energy, we have provided the cen- 
trality dependence of initial conditions at LHC en- 
ergy in Table III and predict the suppression of T 
yields at LHC energy, which is yet to be verified 
with the results available in ALICE experiments at 
LHC energy We have taken three sets of dissocia- 
tion temperatures for the bottomonium states [84| . 
These are obtained by equating a) the lattice free 
energy, b) by subtracting entropy term from the lat- 
tice free energy and c) linear combination of both 
a) and b), with the temperature dependent heavy 
quark effective potential. In Fig. 7, we have plot- 
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FIG. 6: The variation of p T in tegrated survival probability (in the range allowed by invariant p T spectrum of J/ip 
by the Phenix experiment [561 ]') versus number of participants at mid-rapidity. The experimental data (the nuclear- 
modification factor Raa) are shown by the squares with error bars whereas circles and diamonds represent sequential 
melting using the values of To's [5(|] and related parameters from Table II. 



TABLE IV: Formation time (fm), dissociation temper- 
ature To HUl, the speed of sound c 2 a and the screening 
energy density e s (GeV / fm 3 ) calculated in our EoS for 
bottomonium states, respectively. 



State 


TF 


To 


c^(our) 


e s (our) 


T 


0.76 


4.18 


0.322 


496.86 




1.90 


1.47 


0.287 


6.31 


Xb 


2.60 


1.61 


0.294 


9.38 



TABLE VI: Formation time (fm), dissociation temper- 
ature To [13, the speed of sound c 2 s and the screening 
energy density e s (GeV/ fm 3 ), calculated in SQGP EoS 
for bottomonium states,respectively. 



State tf To c 2 (om) e s (our) 

T 0.76 2.90 0.317 111.29 

V 1.90 1.06 0.244 1.47 

Xb 2.60 1.07 0.247 1.58 



TABLE V: Formation time (fm), dissociation tempera- 
ture To [H3], the speed of sound c 2 and the screening 
energy density e s (GeV/ fm 3 ), calculated in our EoS for 
bottomonium states, respectively. 



State 


TF 


To 


c 2 s (our) 


e s (our) 


T 


0.76 


3.40 


0.320 


213.14 




1.90 


1.18 


0.263 


2.44 


Xb 


2.60 


1.22 


0.267 


2.82 



ted the centrality dependence of the p T -integrated 
survival probability for the bottomonium states at 
LHC energy. 



V. CONCLUSIONS 

We revisited the equation of state for strongly in- 
teracting quark-gluon plasma in the framework of 
strongly coupled plasma with appropriate modifica- 
tions to take account of color and flavor degrees of 
freedom and QCD running coupling constant. In 
addition, we incorporate the nonperturbative effects 
in terms of nonzero string tension in the deconfincd 
phase, unlike the Coulomb interactions alone in the 
deconfined phase beyond the critical temperature. 
Our results on thermodynamic observables viz. pres- 
sure, energy density, speed of sound etc. nicely fit 
the results of lattice equation of state with gluon, 
massless and as well massive flavored plasma. Mo- 
tivated by this agreement we apply our equation of 
state to estimate the centrality dependence of J 'ftp 
suppression in an expanding dissipative strongly in- 
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FIG. 7: The variation of p T integrated survival probability versus number of participants for T.The circles and 
diamonds represent sequential melting of rj/s = 1/4tt and rj/s = 0, respectively. The parameter for left, middle and 
right figures given in the Table IV, V and VI, respectively. 



teracting QGP produced in relativistic heavy-ion 
collisions. We have found a complete agreement with 
the PHENIX experimental results on J/^f suppres- 



sion at RHIC energy. Moreover we predicted the 
same for the T suppression which is yet to be veri- 
fied in ALICE experiments at LHC energy. 
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